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ABSTRACT 

We present l.Smilhmeter ALMA Cycle observations of the edge-on debris 
disk around the nearby, ~10Myr-old, M-type star AU Mic. These observations 
obtain 0'.'6 (6 AU) resolution and reveal two distinct emission components: (1) the 
previously known dust belt that extends to a radius of 40 AU, and (2) a newly 
recognized central peak that remains unresolved. The cold dust belt of mass 
~ 1 Mivioon is resolved in the radial direction with a rising emission profile that 
peaks sharply at the location of the outer edge of the "birth ring" of planetesimals 
hypothesized to explain the midplane scattered light gradients. No significant 
asymmetries are discerned in the structure or position of this dust belt. The 
central peak identified in the ALMA image is ~ 6 times brighter than the stellar 
photosphere, which indicates an additional emission process in the inner regions 
of the system. Emission from a stellar corona or activity may contribute, but 
the observations show no signs of temporal variations characteristic of radio- 
wave flares. We suggest that this central component may be dominated by dust 
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emission from an inner planetesimal belt of mass ~ 0.01 Mmooii, consistent with a 
lack of emission shortward of 25 /im and a location <3 AU from the star. Future 
millimeter observations can test this assertion, as an inner dust belt should be 
readily separated from the central star at higher angular resolution. 

Subject headings: circumstellar matter — planet-disk interactions — stars: in- 
dividual (AU Microscopii) — submillimeter: planetary systems 



Introduction 



Debris disks are created by the collisional erosion of planetesimals, the building blocks 
of planetary systems. These collisions continuously generate dust grains with a range of 
sizes that are detected with astronomical measurements from optical to radio wavelengths. 
Resolved observations of nearby debris disks are ins trumental in advan cing our understanding 
of these systems. At a distance of 9.91 ± 0.10 pc f lvan LeeuwenI 120071 ) . the Ml star AU Mic 
hosts one of the closest and best studied debris disks. The detection of submillimeter emission 



(jLiu et al.ll2004j) from this ~10 Myr-old system in the /3 Pic moving group (jZuckerman et al. 



20011) was followed quickly by the discovery of an edge-on disk seen in scattered starlight 
(IKalas et al.ll2004l ). Subsequent work has characterized the scattered light in great detail. 



exploiting its proxiraity to constrain its radial and vertical structure (iLiul l2004l : iKrist et al. 
2OO5I : iMetchev et aPboOsI : ICraham et aPboOTl : IPitzgerald et allboOTh . 



Observations of dust emission at (sub)millimeter wavelengths provide important, com- 
plementary information about debris disk structures. Unlike the small grains probed at op- 
tical and near-infrared wavelengths that react strongly to stellar radiation and wind forces, 
the large grains that dominate the millimeter-wave emission have dynamics more like the 
parent planetesimals. As a result, long- wavelen gth images trace best the location and dis- 
tribution of the larger colliding bo dies (IWvattI 12006.) . and potentially also the signatures 
of planets that interact with them (lErtel et al.ll2012l ). These size-dependent dust dynamics 
manifest beautifully in the edge-on AU Mic disk. Resolved millimeter-wave observations 
show an emission belt within the extended optical disk that pea ks near a radius of 35 AU, 
where the midplane scattered light profile steepens dramatically (jWilner et al.ll2012l ). These 
features are elegantly explained by the presence of a "birth ring" of planetesimals at that 
locati on, where small grains r e leased in a collisional ca scade are launched into an extended 
halo Jstrubbe fc Chianejbood : lAugereau fc BeustI booe l) . 



With the advent of the Atacama Large Millimeter Array (ALMA ), the millimeter e mis- 
sion in nearby debris disks can be imaged in much greater detail (e.g., iBoley et al.ll2012l ). In 
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this Letter, we present new, sub-arcsecond resolution ALMA Cycle observations of AU Mic 
at A = 1.3 mm. The ALMA data provide substantially improved constraints on the locations 
of colliding planetesimals in the AU Mic disk and help shed light on the processes that may 
be shaping the planetesimal distribution. They also reveal a previously unknown, centrally 
located emission feature. 



2. Observations 

AU Mic was observed by ALMA with its Band 6 receivers over four 2 hour- long "schedul- 
ing blocks" (SBs) in 2012 Apr and Jun. Table [T] summarizes the observations. The 16-20 
operational 12-m antennas were arranged to span baseline lengths of 21-402 m (correspond- 
ing to a maximum resolution of ^O'.'G). The correlator was configured to optimize continuum 
sensitivity, processing two polarizations in four 2 GHz-wide basebands, each with 128 spectral 
channels, centered at 226, 228, 242, and 244 GHz. In each SB, we interleaved observations 
of AU Mic (pointing center a = 20M5™09!34, 6 = -31°20'24'.'09, J2000, within 1" of the star 
position at all epochs) with the nearby quasar J2101— 295. 

The data from each SB were calibrated independently within the CASA software package. 
After applying system temperature measurements and phase corrections from the water 
vapor radiometers, the data were flagged and averaged into 6.048 s integrations. A calibration 
of the spectral response of the system was determined from observations of J1924— 292, and 
complex gain variations induced by atmospheric and instrumental effects were corrected using 
observations of J2101— 295. The absolute flux calibration was derived from observations of 
Neptune: a mean calibration was applied to all basebands, with a systematic uncertainty 
of ~10% (see §3.3). To generate an image at the mean frequency, 235 GHz (1.28mm), we 
Fourier inverted the calibrated visibilities with natural weighting and performed a multi- 
frequency synthesis deconvolution with the CLEAN algorithm. The visibilities were further 
reduced by spectrally averaging over the central 112 channels in each baseband and re- 
weighted by the observed scatter. 



3. Results and Analysis 

3.1. Image of 1.3 mm Dust Continuum Emission 

Figure [1] shows an image of the A = 1.3 mm emission from SB-4 (with the most antennas 
and best weather conditions), with synthesized beam 0'.'80 x 0'.'69 (8x7 AU), p. a. 49°, and 
rms of 30 /iJy beam~^. An image constructed from all 4 SBs is consistent but noisier, which 
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Fig. 1. — ALMA image of the 1.3 mm continuum emission from AU Mic. The elhpse in the 
lower left corner represents the O'.'SO x 0'.'69 (8x7 AU) synthesized beam. 



we attribute to systematic calibration issues resulting from the poorer weather conditions of 
the earlier observations. The emission is confined to a narrow band with aspect ratio >10:1, 
with an orientation consistent with the scattered light disk. The emission is not resolved 
in the direction perpendicular to the elongation. There are clear peaks near both extrema 
and in the middle of the structure (detected at all four epochs). The emission is marginally 
brighter at the northwest end than the southeast end, and shows small undulations along its 
length, though none of these variations are significant. We interpret the observed structure 
as a superposition of two components: (1) the nearly edge-on dust belt with limb-brightened 
ansae, and (2) a new, distinct, and compact feature located at the center of the belt. 



3.2. Modeling Formalism 



Building on the phenomenological methodology of IWilner et al.l ( 120111 . |2012[ ) to analyze 
resolved millimeter emission from debris disks, we construct a parametric model to quantify 
the observed properties of the AU Mic emission. We consider two model components: a 
vertically thin, axisymmetric "outer" belt, and an additional source to account for the central 
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peak. The belt component is informed by models of the scattered light that show the disk 
midplane within 50 AU is remarkably straight, <0?5 from edge-on, and thin (FWHM ~0'.'3). 
We assume the belt is viewed at an inclination of 89?5. The belt is characterized as an annulus 
with (unprojected) radial intensity Iu{r) oc for ri < r < ro, with a normalization defined 
by Fbeit = J If dfl, a center determined by offsets (relative to the pointing center) {Aa, A6}, 
and an orientation described by a position angle (PA). We treat the central component as 
a circular Gaussian with mean Arcen, variance a"^^^ (half width at half maximum TZcen = 
v^2 ln2 o"cen); ^ud flux density -Fcen- The mean Arcen is defined as a radial shift from the belt 
center in the plane of the belt. We also include power-law spectral scalings between the 4 
basebands for each component, denoted ctbeit and acen, where F^, oc z/". 

For a given parameter set, we compute four synthetic visibility sets sampled at the same 
spatial frequencies observed by ALMA, corresponding to the spectrally averaged basebands 
(at 226, 228, 242, and 244 GHz). By fitting the visibility data directly, we are not sensitive 
to the non-linear effects of deconvolution, and take advantage of the full range of available 
spatial frequencies. The fit quality is quantified by a likelihood metric, C, determined from 
the values summed over the real and imaginary components at all spatial frequencies 
(ln£ = — x^/2). A Monte Carlo Markov Chain (MCMC) approach was utilized to char- 
acterize the multi-dimensional parameter space of this model and determine the posterior 
probability distributi on functions for each parameter. We used the afiine-invariant ensemble 



sampler proposed by iGoodman fc Weard (|2010|), in a locally - modifi ed version of the paral- 



lelized implementation described by [Foreman- Mackey et al.l (|2012[ ). to compute likelihood 



values for ~10^ MCMC trials. Uniform priors were assumed for all parameters, with bounds 
imposed to ensure that the model was well-defined: {-Fbdt, -F'ccn, o'con) — 0; and < < Tq. 



3.3. Results of Model Fits 

The best-fit parameter values and their 68% uncertainties determined from the marginal- 
ized posterior probability distributions are listed in Table |21 The data and best-fit model 
are compared in the image plane in Figure |2l there are no significant residuals. The best-fit 
model has a reduced = 1-37 (905,920 independent datapoints, 12 free parameters). The 
modeling procedure was performed on each SB individually and the full dataset (all 4 SBs 
together). The results were entirely consistent, although the parameter uncertainties were 
notably smaller from the superior SB-4 dataset alone, and we focus on those results. 

Most parameters are determined with high precision. We find go od agreement of th e 



outer belt parameters {-Fbeit, r^, ro} with the less well-constrained fits of IWilner et al.l (120121) 



and on the disk PA from measurements of scattered starlight (e.g., iKrist et al.ll2005l ). We 
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Fig. 2. — (left) The observed 1.3 mm emission from AU Mic, (center) the best-fit model 
(see §3.3p . and (right) the imaged residuals. Contours are drawn at 4cr (120yuJy beam^^) 
intervals. 



measure a flat spectrum for the outer belt (ttbeit ~ 0) across the 4 basebands, which corre- 
sponds to the difference between the spectral s lopes of AU Mic an d Neptune (ctNeptune ~ 2.1), 
consistent with data from 350 /im to 1.3 mm (jWilner et al.ll2012l ). 

The central emission peak is detected with high confidence at Fcen = 320 /iJy (>10(T 
brighter than the outer belt at that location). It is unresolved, with TZcen < 3.0 AU (3 a), 
and positionally coincident with the outer belt center: Arcen < 1.9 AU (3 a). Regarding the 
outer belt, the most notable result is that the models strongly favor rising emission profiles 
with large, positive gradients: x ~ 2.3 ± 0.3. Models with the standard assumption of x < 
produce significant residuals, under-predicting the intensities at ±1-2" from the belt center. 
Because of the steep increase in the emission profile, there is only a weak constraint on the 
inner edge of the outer belt. The best-fit rj deviates from at the ~2(T level: the 3 a limit 
is ri < 21 AU. 



4. Discussion 



We have presented new, sub-arcsecond resolution ALMA observations of 1.3 mm emis- 
sion from the AU Mic debris disk and analyzed the data with a simple parametric model. 
This emission is resolved into two distinct components: (1) an edge-on outer belt with an 
emission profile that rises with radius out to 40 AU, and (2) an unresolved peak at the 
center of the outer belt. This distribution is more complex than the single, narrow ring often 
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assumed for debris disks. Howe ver, it has son ae similarities t o other nearby resolved systems, 
like e Eri ( Backman et al.ll2009l ) or HR 8799 JSu et alibopgl ). that show an inner component 
inferred from excess infrared emission, separate from an extended and colder outer belt. 



4.1. The Central Emission Peak 



The stellar photosp here is much fainter th an the central peak noted in Figure [H A 
NextGen stellar m odel faauschildt et all Il999h wit h T.ff = 3720 K, = 0.11 L©, and 



= 0.6 M0 (e.g., iMetchev et all 1200,4 IChen et al.ll2005h that matches the AU Mic pho- 
tometry from 0.4-25 fim contributes only F^, = 52 /iJy at 1.3 mm, ~6x fainter than observed. 
However, AU Mic is an active star that exhibits radio- wave bursts. In quiescence, observa- 
tions find <120 /iJy at 3.6 cm ( IWhite et al.lll994j ). and the contribution at 1.3 mm from hot 
coronal plasma se en in X-rays is un likely to be significant (though better spectral constraints 
are d esirable, see iLeto et al.ll2000l ). Flares are detected from AU Mic at ~200-1200 /iJy at 
6 cm jBower et al.ll2009l ). but this non-thermal emission is much weaker at 1.3 mm. While 
the unknown variability makes any extrapolation to 1.3 mm problematic, the temporal prop- 
erties of the ALMA emission p rovide additional in formation. Radio-wave flares have fast 
decay times, of order an hour (IKundu et al.lll987l ): but, the mm-wave peak persists at a 
consistent intensity in all four ALMA observations, within uncertainties that are typically 
2 — 3x larger than for SB-4, spanning timescales from 1 hour (within SB-4) to 2 months 
(SB-1 to SB-4). Unfortunately, the spectral index (ocen) constraints are not good enough to 
be diagnostic. We suspect that stellar emission is too weak and too ephemeral to be respon- 
sible for the 1.3 mm peak, but the available data does not allow for a firm determination of 
its contribution. 

Alternatively, the central emission peak could be produced by dust in a distinct (un- 
resolved) planetesimal belt located close to the star. In §3.3, we constrained the extent of 
this peak to TZcen < 3 AU (3cr), inside the inner working arigle (0?8 ^ 8AU) of all previous 



high resolution imaging of scattered light (IKrist et al.l l2005l : iFitzgerald et al.l 120071 ) . Rough 
models of the spectral energy distribution (SED) from the ALMA central peak can help 
assess the feasibility that it originates in an inner dust belt. In this context, the m ost salient 



feature of the AU Mic SED is the absence of emission excess at A < 25 /xm (e.g., iLiu et al. 



2004J : IChen et al.ll2005l ). We assume the central peak represents the combined emission from 
the star and dust, such that -Fdust = -^cen — -F* ~ 0.25 mJy at 1.3 mm. Optically thin dust 
emission at a temperature, T, has Fdust ~ i^uBiy(T)Mdnst/cP, where is the opacity spec- 
trum, the Planck function, Mdust the mass, and d = 9.91 pc. For a given dust population 
characterized by Kiy, we computed the maximum T (and minimum Mdust) consistent with 
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both the observe d miUimeter flux density and the infrared SED. We calculated various K,y 
for dust with the IWeingartner &: Draind (120011 ) "astrosilicate" compositi on and a power-law 
size d istribution n{a) oc a~^-^ between a^i^ = 0.2 /im (the blow-out size; IStrubbe &: Chiang 
20061 ) and Omax values from 1 fim to 1 cm. For Omax < 100 /xm, models of the central peak 



over-predict the observed 60-70 /zm emission if T > 35 K. However, larger grains with 
Omax > 1 mm at temperatures up to T 75 K can be accommodated without producing 
an excess at A < 25 fim. These maximum T values are comparable to the expected dust 
temperatures a few AU from the star, compatible with the emission size constraints (T^cen)- 
The corresponding minimum Mdust is ~9 x 10^^ g, about 1% of the lunar mass. These calcu- 
lations show that the central emission peak is consistent with a cool dust belt located <3 AU 
from the central star, with a total mass comparable to the asteroid belt in our Solar System. 
If this interpretation is correct, then ALMA observations at higher resolution can determine 
its properties. Interestingly, the temperature of this putative inner belt is c older than the 



y 190 K found to apply systematically to inner belts around F5-K0 stars by [Morales et al. 
J201lh . 



4.2. The Outer Dust Belt 



Our modeling of the ALMA data locates the far edge of the outer emission belt with 
high precision, = 40 AU, which matches closely the outer edge of the hypothesized "birth 
ring" of colliding planetesimals. This analysis does not define the shape of the edge below the 
~ 6AU resolution limit, but the truncation is re miniscent of the outer edge of the classical 
Kuiper Belt (47 ± 1 AU; iTrujillo fc Brown! 1200 ll ). T he origins of su c h sharp edges rem ain 
unclear: they could be from dynamical interactions (jida et al.ll2000l : iBoley et al.l 120121 ). or 
they may simply represent the initial conditions, where planetesimal formation was efficient 
and successful in the primordial disk. Adopting the opacity used in §4.1 {ki^ = 2.7cm^g~^), 
and assuming T ^ 25 K ( for 35-45 AU), t he dust mass of this outer belt is 7x 10^^ g (consistent 
with previous estimates; iLiu et al.ll2004l ). ~100x more massive than the hypothesized inner 
belt; the Kuiper Belt and asteroid belt have a similar mass ratio. 

The mm-wave emission morphologies of cold belts of dusty debris reflect the dynamical 
processes that shape the underlying planetesimal distributions. For AU Mic, our modeling 
suggests that its outer emission belt can be described by an increasing emission proflle with 
a positive radial power-law index x ~ 2.3 ±0.3. If we assume the emitting dust is in radiative 
equilibrium with a temperature proflle T oc r~°'^, this implies a rising surface density proflle, 
S oc r^'^, strongly peaked near 40 AU. A broad parent body ring with constant surface 
density would produce a radial intensity proflle with x ~ —0.5, a value ruled out with 
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high confidence ( > 5a). A rising behavior is predicted for "self-stirred" disks with ongoing 
planet formation (IKenyon &: Bromleyll2002l ): in particular, the models of Kennedy fc Wyatt 
(120 lOl ) suggest S oc r^/^. However, the timescale required to assemble Pluto-sized bodies at 
~40 AU to initiate a c oUisional cascade around a low- mass star like AU Mic is much longer 
than its ~10Myr age (IKenyon fc Bromleyll2008l ). Moreover, this scenario does not naturally 
accommodate the presence of a separate, interior planetesimal belt. Of course, the still 
modest resolution of the data is compatible with more complex scenarios, such as multiple 
closely-spaced belts of different brightnesses that mimic a smooth gradient. Scattered light 
observations of the AU Mic disk show asymmetries on both large and small scales, with 
several peaks and depressions projected against the broad ansae in Fig ure [H at radii beyond 
the millimeter undulations (features A-E; see Fitzgerald et al.l 120071 ). With such a steep 
emission gradient in this outer belt, the data do not strongly constrain its width, or the 
location of its inner edge. Our modeling indicates substantial emission from mm-sized grains 
interior to 40 AU, in the ^20- 40 AU zone inferred to be highly depleted of /im-sized grains 
from polarized scattered light (IGraham et al.l 120071 ). 



The ALMA data show no clear evidence for asymmetries or substructure that would 
signal planet-disk interactions. The hints of modulating millimeter brightness along the belt 
in Figure [1] are insignificant in the residuals from subtracting a symmetric parametric model 
(see Figure [2]). This rules out substructure brighter than 90 /xJy beam~^ (3 a), corresponding 
to dust clumps > 1% of the lunar mass (for the dust properties adopted above). Those 
limits argue against over-densities of dust-producing planetesimals trapped in mean motion 
resonances (IKuchner fc Holmanll2003l ). as might arise from the outward migration of planets 
Jwvattlbooah . Given the young age of the system, the broad and smooth character of 
the outer belt in the AU Mic disk may resemble the Kuiper Belt prior to the epoch of 
Neptune's migration (jMalhotralll995l ) . It is interesting that none of the cl aims of millimeter 



emission clurnps in debris disks have survived scrutiny at higher sensitivity (jPietu et al.ll2011 



Hughes et al. 


2011 


2012) 


(Kuchner & StarkI 


2010|). 



20121) . It may be that any such features are effectively erased by collisions 



and c entral peak, as mig ht result from the secular perturbations of a planet in an eccentric 
orbit (IWyatt et al.lll999l ). The limit on the displacement, Arccn < 1.9 AU (3 a), corresponds 



approximately to a limit on ae, where a is the semi-major axis and e is the eccentricity. 
This limit can still accommodate a wide-orbit planet with modest eccentricity, similar to 
Uranus. Such a planet could be responsi ble for stirring the disk to 40 AU in ~10 Myr (e.g. 
for a = 30 AU and e = 0.05, see eqn. 15 of lMustill &: Wvattll2009l). Limits from high c ontrast 



direct imaging admit Saturn-mass planets at these separations (jPelorme et al.ll2012l ). 
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4.3. Concluding Remarks 

The basic architecture of the AU Mic debris disk appears remarkably similar to the 
Solar System, with a potential analog to the asteroid belt at a few AU, and a colder, more 
massive, and apparently truncated counterpart of the Kuiper Belt extending to 40 AU. Future 
observations are needed to determine if stellar processes could be responsible for emission 
attributed to the asteroid belt, and to determine if the Solar System analogy extends to 
include a planetary system like our own. 
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NAOJ. The National Radio Astronomy Observatory is a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 
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Table 1. ALMA Cycle Observations of AU Mic 



ID 


Date (UT) 






Antennas 


PWV (mm) 


SB-1 


2012 Apr 23 07:30 - 


09 


26 


17 


1.7 


SB-2 


2012 Apr 23 09:39 - 


11 


03 


16 


1.7 


SB-3 


2012 Apr 24 09:09 - 


11 


19 


18 


3.0 


SB-4 


2012 Jun 16 05:48 - 


08 


02 


20 


0.7 



Table 2. Model Parameters 



Parameter 


Description 


Best-Fit 


68% Confidence Interval 


-fbelt 


belt flux density (mJy) 


7.14 


+0.12, -0.25 


X 


belt radial power law index 


2.32 


+0.21, -0.31 


Ti 


belt inner radius (AU) 


8.8 


+11.0, -1.0 


To 


belt outer radius (AU) 


40.3 


+0.4, -0.4 


PA 


belt position angle (°) 


128.41 


+0.12, -0.13 


Qibelt 


belt spectral index 


-0.15 


+0.40, -0.58 


^ cen 


Gaussian flux density (mJy) 


0.32 


+0.06, -0.06 


Arcen 


Gaussian offset (AU) 


0.71 


+0.35, -0.51 


2 


Gaussian variance (AU^) 


<5.9 


(3 a limit) 




Gaussian spectral index 


-0.35 


+2.1, -4.5 


Aa 


R.A. offset of belt center (") 


0.61 


+0.02, -0.02 


A5 


Dec. offset of belt center (") 


-0.03 


+0.02, -0.02 
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